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In this communication, the mass declined for the decay channel, e+e 



ZH 



e e + X, as measured by the ILD detector was studied. The Higgs mass is assumed 
to be 120 GeV and the center of mass energy is 250 GeV. For an integrated luminosity 
of 250 fb _1 , the accuracy of the reconstruction and the good knowledge of the initial 
state allow for the measurement of the Higgs boson mass with a precision of about 
100 MeV. 



1 Introduction 

One of the major goals of the ILC (" International Large Collider" ) is to observe the Higgs bo- 
son and describe its properties. In the Standard Model [I] the Higgs particle plays a key role 
in explaining the origin of the masses of the 
elementary particles. If existing, a Higgs 
boson with a mass of 120 GeV is favoured by 
recent analyses of electo-weak data [3] . The 
Higgs-strahlung process (e + e~ — s- Z* —> 
ZH) is the major Higgs production mech- 
anism at ILC, see Figure [TJ In particular, 
in this study, we will consider the channel 
e+e- -> ZH -> e+e- +X. 

By detecting the e+e~ lepton pair pro- 
duced by the well-known Z boson decay, 
one can measure the Higgs mass (using the 
mass recoiling of the Z) independently of 
the Higgs decay. The most favourable case 
is obtained for the center of mass energy 

slightly larger than the Z boson mass plus the Higgs boson mass (i.e. 




Figure 1: Higgs boson production via Higgs- 
strahlung process. 



211.2 GeV) 0. 



2 Event generation and detector simulation 

The analysis was performed with the International Large detector (ILD) data sample ful- 
filled for ILD Letter of Intent [4] . The events are fully simulated and reconstructed with the 

* Research conducted in the scope of a LIA known as "International Laboratory for Collider Physics - 
ILCP". The moroccan contribution to this work is also supported by the High Energy Physics Network 
(RUPHE). 
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ILD_00 detector model. The Higgs mass is assumed to be 120 GeV, and the center of mass 
energy is 250 GeV, with an integrated luminosity of 250 fb -1 . The e + polarization is taken 
to be equal to 30% and the e~ polarization is 80%. 

Table Q] gives the signal (ZH— >eeX) cross 



Cross section 
(fb) 


c + c beam polarization 
mode (30%, 80%) 


Process 


(+,-) 


(-,+) 


ZH^eeX 


12.55 


8.43 


ee (Bhabha) 


17.30 10 B 


17.28 10 s 


4f->ccfT 


4897 


3793 



section as a function of the beam polar- 
ization (e.g. (+,-) mode means that the 
positron beam polarization is +30% and the 
electron beam polarization is -80%). The 
polarization modes (-,-) and (+,+) lead to 
smaller signal cross sections: 7.65 and 5.78 
fb respectively for similar levels of back- 



Table 1: ZH^ccX and background cross sec- ground. The Bhabha scattering (ee) and the 
tions as a function of the beam polarization. Standard Model events with 4 fermions, in- 
cluding e + c~ contribute to the background. 
The cross section for Bhabha scattering is 6 orders of magnitude larger than the signal cross 
section. 



3 Event reconstruction and background rejection 



The identification of the Z boson is obtained by selecting the e + e pair which gives the best 
mass for the Z boson. In our case: y^zreconstructed = M^i 10 GeV. In addition, central 



leptons with opposite charge are required, i.e. 
measurement is optimal) 



(+,-) 


(fb" 1 ) 




N rcco » at 


^expected 
for 250 fb -1 


ccX 


10 000 


125 500 


49 799 


1 245 




0.5123 


S S66 734 


48 201 


23 10 6 



Table 2: Number of simulated, reconstructed 
and expected events with their associated lu- 
minosity for the (+, — ) beam polarization 
mode. 



|cos(6>.;)| < 0.9 (in that case the momentum 



The Higgs mass is then calculated from 
the well known formula: 

M 2 H = s + M| - 2M Zy /s, 
where Mz is the reconstructed mass 
of the Z boson and yfs the total en- 
ergy in the mass center of the colli- 
sion. 

The expected number of events after selec- 
tion is given in Table [5] Despite a wor- 
thy gain, the Bhabha effect remains strongly 
dominant. 



4 The Bhabha scattering 

To reject the Bhabha events, a cut on the number of reconstructed objects (J^objects > 21) 
is applied, see Figure [5] This cut removes the first maximum in the signal distribution 
which probably corresponds to events where the Higgs boson decays into a t + t~ pair. 
Unfortunately, the measurement of the cut efficiency is limited by the available Monte Carlo 
event number, see Table Nevertheless, this cut looks very efficient because of the shape 
of the Bhabha distribution. Therefore, in the following parts the Bhabha effect will be 
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neglected but the cut on the number of reconstructed objects will be maintained. 



Signal 
Background 



(+,-) 




^objects > 21 


^objects > 21 




(fb- 1 ) 




for 250 fb" 1 


eeX 


10 000 


55 847 


1396 


ee 


0.5123 


1 


487 




Table 3: Number of events, with their as- 
sociated luminosity, after the cut on the 
number of reconstructed objects for the 
(+, — ) beam polarization mode. 



Figure 2: Normalized distributions (surfaces 
equal to 1) of the reconstructed objects from 
the Particle Flow Algorithm (PFA) [5] for the 
ZH signal (in blue) and the ee Bhabha scat- 
tering (in red). 



5 Event selection (cuts on kinematic variables) 

The following cuts were applied to reject the 4 fermion background: 

1. the di-lepton transverse momemtum calculated from the vectorial sum of the two 
leptons: 18 GeV < Pt dj/epto „ < 68 GeV, 

2. the acollincarity (i.e. angle between the 2 leptons, acol = cos _1 (Pi.P2/|Pi|.|P2|) ): 
0.4 < acollinearity < 1.35, 

3. the energy of the Z boson: 20 GeV < E z < 115 GeV, 

4. the Higgs mass interval: 115 GeV < M H < 165 GeV. 

Different cuts are shown in Figure[3] Reconstituted recoil mass distributions before and after 
the selection are presented in Figure [4] 
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E z (GeV) ~ m_H(GeV) 



Figure 3: Event selection cuts (the distribution surfaces are normalized to 1). 




Figure 4: Reconstruted recoil mass distributions before and after the final cuts. 
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6 Reconstruction of the Higgs mass 

The Higgs recoil mass distribution is fitted with a GPET function ("Gaussian Peak Expo- 
nential Tail" ) defined in the following way: 



f(x\ a, n, x,a) = N 



(x-x) 2 

e : for ^—^ < a 

J (J — 

P . e ~ {J ^- + (f -/?) -e^r ■ e -(*-*)-§ : for ^ > a 



where /3 € (0,1)- This function and its first derivative are continuous when (x — x)/a = a. 

When the cuts are applied to reject the background, one observes a small distorsion in 
the signal distribution in the region between 121 and 123 GeV. This can be removed by 
adding an additional Gaussian in the same region, as given by: 

/ (x; a, n, x, a) = f(x; a, n, x , a) + NiG{x;x\,ui) 

The parameters N\, x\ and o\ are determined in the fit {x\ is around 122 GeV and o\ is 
between 1 to 3 GeV). The corrected GPET function will be used in the following studies. 

Figure [5] shows the pure signal fitted by the corrected function for the (+,-) polarization 
mode. 

The background fits to a polynomial sum of degree 6, see Figure [SI 




Figure 5: Fit of the Higgs recoil mass dis- 
tribution with the corrected GPET function 
for the (+,-) beam polarization mode. 



Figure 6: Fit of the background distribution 
with a polynomial sum for the (+,-) beam 
polarization mode. 
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7 Measurement of the Higgs mass with signal and background 



Figure [7] shows the Higgs recoil mass distributions (signal + background) for different beam 
polarization modes. Some of the fit parameters are given in Table [4] The accuracy on the 
Higgs mass measurement is around 1 00 McV. 




115 120 125 130 135 140 145 150 115 120 125 130 135 140 145 150 

Hecoii(GeV) M recoil (GeV) 

Figure 7: Fits of the Higgs recoil mass distributions (signal + background) for different 
beam polarization modes. 





e + e beam polarization mode (30%, 80%) 


(+r) 


(-,+) 


Only signal 


Signal + Bkgrd 


Only signal 


Signal + Bkgrd 


M H (GcV) 


120.486 ± 0.073 


120.368 ± 0.100 


120.507 ± 0.085 


120.445 ± 0.110 


a (GeV) 


0.638 ± 0.051 


0.575 ± 0.083 


0.654 ± 0.062 


0.592 ± 0.100 



Table 4: Fit parameter results of the Higgs recoil mass for different polarization beam modes. 
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8 Conclusion and outlook 



The measurement of the Higgs mass recoil in the channel: e + e~ — > Z H — > e + e~ + X 
illustrates the ILC potential for accurate measurements. The effect of the background is to 
deteriorate the accuracy on the Higgs boson mass. For Mh = 120 GeV and a luminosity of 
250 fb _1 , it was found: 

Mfl ec = 120.368 ± 0.100 GeV, for the beam polarization mode (+,-) and 
M Rec = 120.445 ± 0.110 GeV, for the beam polarization mode (-,+). 

The Higgs recoil mass distribution is asymmetric. This could be, for instance, partially 
corrected by taking into account the bremsstrahlung photons emmitted in the ILD detector 
and by using the electromagnetic calorimeter and other sub-detectors. This is the next step 
in the analysis. 

For more information please refer to [6l [2] . 
Acknowledgments 

The authors would like to thank Mr. H. Li for his help and his numerous advice. 
References 

[1] S. Weinberg, Phys. Rev. Lett. 19 (1967) 1264. 

A. Salam, Elementary particle theory, edited by N. Svastholm (Almquist et Wiksell, Stockholm, 1968). 
S.L. Glashow, Nucl. Phys. 22 (1961) 579. 

[2] H. Li's thesis: "Higgs Recoil Mass and Cross-Section Analysis at ILC and ..." 

http: / /tel.archivesouvertes.fr/index.php?halsid=3t3g7n01u9vli57j7u7b5fd763 feview_this_doc=tel- 
00430432&version= 1 

[3] The LEP Elcctroweak Working Group, arXiv:0811.468 2][hep-ex] (November 2008) and up dates for 2009 
summer conferences, see http://lepcwwg.wcb.ccrn.ch/LEPEWWG/plots/summcr2009/ 

[4] The International Large Detector, Letter of Intent, by the ILD Concept Group, March 2009 
|http:/ / www.ilcild.o rg/documents/ild-lcttcr- of-intentl 

[5] M.A. Thomson, Particle Flow Calorimetry at the ILC, Pramana 69 (2007) 1101, |arXiv:physics/0607261| 
[physics, ins-dct] . 

[6] "HZ Recoil Mass and Cross Section Analysis in ILD", H. Li, R. Poschl and F. Richard, LAL 09-121, 
LC-PHSM-2009-006. 



LCWS/ILC2010 



